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ABSTRACT 


A  plasma  simulation  technique  has  been  developed  which  can  be  used 
to  study  the  effects  of  homogeneous  or  inhomogeneous  plasma  sheath  on  the 
radiation  pattern  deterioration  and  input  impedance  of  microwave  antennas. 

A  tank  has  been  designed  and  constructed  for  use  in  the  simulation  technique 
which  can  reproduce  by  means  of  real  dielectric  materials  the  dielectric 
constant  encountered  in  plasma  covered  antenna  research. 

The  radiation  patterns  and  the  input  impedances  of  an  annular  slot 
and  a  thin  and  long  rectangular  slot  have  been  successfully  measured  in  the 
presence  of  a  simulated  loss-less,  homogeneous  and  isotropic  plasma  layer  of 
varied  thickness.  Comparison  with  the  available  theoretical  data  indicates 
a  generally  good  agreement,  although  some  differences  exist.  In  the  case  of 
the  radiation  patterns  these  differences  are  attributed  to  the  finite  distance 
between  the  radiator  and  the  receiving  antenna  on  one  hand  and  the  inherent 
inaccuracy  of  the  saddle  point  method  of  integration  in  certain  region  3  on 
the  other  hand. 
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1.  Introduction 


Recent  aerospace  achievements  have  motivated  a  great  number  of 
investigations  into  wave  propagation  through  plasma  sheaths.  Most  of  these 
investigations  are  of  theoretical  nature  and  they  use  an  infinite,  isotropic, 
homogeneous,  and  uniform  plasma  slab  as  a  model  for  reentry  environment  of  a 
radiating  system.  In  order  to  experimentally  verify  the  theoretical  work  and 
to  obtain  the  results  for  more  realistic  radiating  system  and  plasma  sheath 
configurations,  a  laboratory  simulation  of  a  plasma  sheath  is  very  desirable. 
The  most  recent  simulation  attempts  make  use  of  artificial  dielectrics,  especi¬ 
ally  rodded  media  [kotman,  1962;  Golden  and  Smith,  1964;  Golden,  1964]].  These 
attempts  are  based  on  the  fact  that  rodded  media  have  refractive  indices  less 
than  unity  under  certain  conditions  [Brown,  1953] . 

The  results  of  these  experiments  generally  match  the  theory.  Certain 
important  points,  however,  were  not  found  to  coincide  with  the  theory.  For 
example,  the  sharp  peaks  at  the  critical  angle,  sin  6C  =  n,  were  not  found  as 
expected  when  an  E  plane  slot  antenna  was  used  [Golden,  1964]  .  Such  discre¬ 
pancies  are  possibly  due  to  "higher-order  reflections"  from  the  rodded  media 
which  result  with  increasing  oblique  incidence  [Brown,  1953]  . 

It  is  the  purpose  of  this  report  to  introduce  another  method  of 
simulating  the  plasma  sheath.  This  method  will  be  seen  to  provide  an  exact 
determination  of  the  simulated  plasma  parameters  and  a  laboratory  geometry  for 
which  it  is  a  simple  matter  to  accurately  adjust  and  vary  the  slab  thickness. 

The  simulation  method,  to  be  described  in  detail  in  a  later  section, 
permitted  the  recording  of  the  field  patterns  in  the  region  outside  of  the 
plasma  sheath.  The  radiators  used  were  an  annular  slot  and  a  rectangular  slot 
in  a  ground  plane.  Field  strength  and  admittance  measurements  were  made  for 
each  antenna  with  various  slab  thicknesses.  The  comparison  of  the  experimental 
results  with  theory  determines  the  feasibility  of  this  method  as  a  valid  plasma 
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sheath  simulation  technique.  This  feasibility  has  been  established  for  a 
lossless  plasma  and  in  the  case  of  thicker  plasma  slabs  it  is  limited  only 
by  the  size  of  the  simulation  tank  used. 

2.  Theoretical  Background 
2.1  Annular  Slot  Antenna 

An  annular  slot  cut  in  an  infinitely  conducting  plane  and  covered 
by  a  thin  layer  of  lossless  isotropic  and  homogeneous  plasma  will  be  consider¬ 
ed  first.  The  geometry  of  the  problem  is  shown  in  Figure  1. 

The  slot  is  assumed  to  be  of  narrow  width  with  a  mean  radius  "b". 

A  constant  voltage  "V"  is  assumed  to  be  present  across  the  slot.  The  hori¬ 
zontal  plane  z  »  a  coincides  with  the  interface  of  the  homogeneous  plasma 
layer  and  the  surrounding  medium.  For  convenience,  reference  will  be  given 
to  the  surrounding  medium  (1)  and  plasma  medium  (2)  with  relative  dielectric 
constants  and  ^  y  respectively. 

It  is  clearly  seen  from  the  symmetry  of  the  problem  that  there  is 
no  variation  in  the  ^-coordinate,  i.e.,  ^  «  0  ,  and  only  the  component 
of  the  magnetic  field  exists.  Furthermore,  it  is  noted  that  the  electric 
field  at  the  aperture  is 


Ep  =  -  V5(p  -  b)8(z) 


(1) 


Assuming  the  e”*wt  time-dependence  and  suppressing  it  throughout,  the  electric 
field  components  are  given  by 


1  cM* 

2  s  1  ■  ^ — 

p  iwe0e  oz 


and 


s* =  -  p  h  (pH*> 


(2a) 


(2b) 


where  satisfies  the  equation 


5p 


i 

Lp  c>p 


(pH^) 


d2H, 

&z2 


*  +  k2H$ 


-V5(p-b)5(z) 


0  <  z  <  a 


0 


z  >  a  . 


(3) 


The  solution  of  (3)  has  been  found  to  be  [Tyras,  1962] 


3 


H*1 


iVbtoe0nz 

~  i 


e->,l(z-a)j1(Xb)Hi1\xp)  ‘  XdX 
[Y2  sinh(y2a)  +  n^  cosh(Y2a>] 


(4a) 


iVbwe0n2 
H42  =  2 


[Y2COsh  Y2(a"z)  +  n^Yisinh  Y2(a“2)] 
Y2[y2  sinh(Y2a)  +  n2Y1  cosh(Y2a)] 


JiCXb)^^  (Xp)XdX 
(4b) 


where 

y±  =  \/x2  -  k2  (4c) 

«  2n/Ki  ,  n 2  *  62^1  (4d) 

Changing  to  spherical  coordinates  in  the  configuration  space  and  taking 
a  -  k-^sinp  in  the  transform  space,  (4a)  is  changed  into  a  contour  inte¬ 
gral  in  the  complex  p-plane.  The  evaluation  of  the  contour  integral  by 
the  saddle  point  method  gives  the  far-zone  field 


Y^  ^  n2  b  J^(k^b  sin0)  cos0 

i/n2  -  sin20  sin(k^a  V  n2  -  sin20  )  +  in2  cos0  cos(k]^a  v^n2  -  sin20  ) 


•  ej&jlr  -  n/4]  (5) 

r 

where  Y^  =  \j  €0£i/p0  is  the  admittance  of  the  surrounding  medium.  The  far- 
zone  electric  field  has  the  single  component 

E0i  “  zi  Hn  (6) 

where  Z]Y^  =  1  .  It  is  clear  that  (5)  is  valid  for  real  as  well  as  complex 
values  of  the  plasma  index  of  refraction,  n  , 
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The  expression  for  the  far-zone  electric  field  has  been  programmed 
for  the  computer  [Tyras,  1962]  for  different  combinations  of  frequency  and 
plasma  thickness  and  the  results  are  shown  in  Figure  2.  The  apparent  feature 
that  is  common  to  all  of  these  radiation  patterns  is  the  sharp  attenuation 
for  all  angles  larger  than  the  critical  angle.  This  fact  is  well  explained 
from  the  consideration  of  geometric  optics.  Moreover,  it  is  also  evident  that 
with  the  increasing  electrical  thickness  of  the  plasma  layer,  progressively 
more  ripples  appear  in  the  radiation  pattern. 

Equation  (5)  gives  the  complete  solution  of  the  fields  in  the  radia¬ 
tion  zone.  However,  “leaky  wave"  poles,  poles  arising  from  the  singularities 
of  the  contour  integral  and  located  in  the  improper  Riemann  sheet,  may  signifi¬ 
cantly  affect  the  near  and  intermediate  zones.  These  "leaky  waves"  exhibit  a 
radial  attenuation  and  their  effect  on  the  radiation  field  diminishes  even 
when  strongly  excited.  An  extensive  treatment  of  complex  poles  is  given  by 
Tamir  and  Oliner  [l962], 

2.,  2  Rectangular  Slot  Antenna 

The  second  kind  of  antenna  considered  is  a  thin  rectangular  slot 
in  an  infinite  ground  plane  clad  with  plasma.  The  geometry  of  the  problem  is 
shown  in  Figure  1.  Assuming  that  in  the  aperture  there  is  a  single  component 
of  the  magnetic  field,  ,  then  it  follows  that  the  electric  fields  will  be 
given  by 


Ey 


l 


5T 


(7a) 


and 


-1 

u>e0e  §y~ 


(7b) 


The  solution  can  be  formulated  in  a  straightforward  manner  using 
Fourier  transforms  {knop  and  Cohn,  1964]  .  One  obtains 

W6on  pr\ 

%i  *  ~4^r  JJ  ^y(^n)d^n 

ap. 


CO 

,  IT  ej [gifr-e)-«a(y-»i)-Ki(»-«)1  dc(lda2 

n2s^cos(s2a)  -  is2sin(s2a) 


(8) 


for  the  air  region  and 


weon?'  rr 

nx2  =  “^5“  Jl  Fy <6 

ap. 


00 

If 


in2sxsin(s2(z-a))+S2COs(s2(z-a)) 

o  - -  -  — - 

S2(n^s^cos(s2a)-is2sin(s2a)) 


ei[«L  x*’0+°!2(y“Tl)]  dOfjdC^ 


(9) 


for  the  plasma  region,  where 

st  =  \l kj  -  (c^2  +  C^2) 


(10) 


and  Fy(|,n)  is  the  y-th  component  of  the  electric  field  in  the  aperture  which 


is  assumed  to  be  known. 
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For  the  particular  case  at  hand  when  the  slot  is  narrow  and  very 
long  one  can  assume  the  excitation  of  the  form 

*y(^)  «  -V8(n)  (10) 

where  V  is  the  voltage  across  the  slot.  Substituting  (10)  into  (8)  and 
using  the  saddle  point  method  of  integration,  one  obtains  for  the  radiation 


field 

we  n2V  i(kop-n/4) 

Hxo  *  ~f=~  ■  G(8)  S~F= - 

V  2jt  V  k0p 

where  G(0)  is  the  form-factor  given  by 


(11a) 


G(0) 


cos  0 


n^cos  9  cos  (koa\/n^-sin^0)  -  iv4i^-sin^0  sin  (kQal/ n^  -  sin20) 

(lib) 


The  expressions  in  (11)  have  been  derived  earlier  by  Newstein  and 
Lurye  [l956]  ,  Omura  [l962]  and  Tamir  and  Oliner  [l962].  Omura' s  report 
contains  a  considerable  number  of  numerical  results  to  which  we  shall  later 
have  the  occasion  to  refer  to  for  comparison  with  our  experimental  results. 
The  work  of  Tamir  and  Oliner  is  devoted  to  interpretation  of  the  results  in 
terms  of  "leaky  waves". 

Both  the  numerical  results  of  Omura  [l962]  and  the  "leaky  wave" 
analysis  of  Tamir  and  Oliner  [l962],  as  well  as  the  simple  consideration 
of  geometric  optics,  predict  that  the  radiation  pattern  will  have  peaks  in 
the  neighborhood  of  the  critical  angles.  Moreover,  as  the  thickness  of  the 
layer  increases,  progressively  more  minor  peaks  appear  in  the  radiation 
pattern  while  the  major  lobes  are  still  found  ir.  the  vicinity  of  the  critical 
angles  [Tamir  and  Oliner,  1962]  . 


3.  The  Experimental  System  and  Results 
3.1  The  Simulation  Technique 

Plasmas  are  characterized  by  a  dielectric  tensor  which  reduces 
to  a  complex  scalar  in  the  absence  of  external  magnetostatic  fields. 
Furthermore,  the  real  part  of  the  plasma  dielectric  constant  is  less  than 
unity,  namely,  ep/ea  <  1  where  ep  is  the  plasma  dielectric  constant 
and  ea  is  the  dielectric  constant  of  free  space.  If,  however,  ep/ea  * 
epS/eas  <  1  where  the  subscripts  nps"  and  "as”  denote  plasma  simulation 
and  air  simulation,  respectively,  it  is  seen  that  a  simulated  plasma 
environment  depends  on  the  ratio  of  the  dielectric  constants  and  not  on 
their  absolute  values.  Thus  an  artificial  plasma  environment  can  be 
created  by  covering  a  radiator  or  a  scatterer  under  investigation  with 
a  medium  having  a  dielectric  constant  less  than  that  of  the  free  space 
simulator. 

A  major  requirement  of  the  free  space  simulator  is  that  it  be 
in  the  liquid  form  to  allow  movement  of  measuring  equipment  through  it. 
Furthermore,  it  should  have  a  low  loss  tangent  to  minimize  signal  attenu¬ 
ation  in  the  medium.  The  plasma  sheath  simulator  can  be  any  foamy 
material  with  er  ~  1  or  simply  air  itself.  Such  a  combination  of  di¬ 
electric  materials  will  simulate  a  plasma  with  0  <  ep  <  1  . 
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With  such  a  simulation  technique,  the  plasma  parameters  can  be 
scaled  and  properly  defined.  Consider  the  annular  slot  as  an  example. 

The  far  zone  electric  field  of  the  annular  slot  covered  with  the  dielectric 
layer  can  be  expressed  from  (5)  as 

bVk  n^  J,  (k  b  sin0)  cos 0  e^^^°r  ”  •  r 

„  O  JL  O 

Eq  *  -  ■  -  — -  — - -  —  ■  — 

sin^S  sin(koaV^n^-sin^0)+in^cos0cos  (koa\/n^-sin^0) 

(12) 

on 

where  n  *  1  -  (Wp/u)^  ,  "a"  is  the  thickness  of  the  dielectric  layer,  and 
"b"  is  the  mean  radius  of  the  slot.  If,  however,  the  air  is  replaced  by 
another  dielectric  medium,  as  in  the  proposed  simulation  technique,  the 
far~zone  field  takes  on  the  form 


b  Vkn2  J.(kb  sin0)  cos 0  e^s1'*/4)  .  r-1 
„  sssl  ss 

Eg  =*  -  .  .  - .  .  —  - -  ...  - - ~  '■■  ■  —  ■  - 

Vn^  -  sin^0sin(k  a  l/n^-sin ^0)+  in^cos0cos(k  a  'I  n2  -sin^0  ) 
s  '  s  s  s  sss  • 

(13) 

where  the  subscript  "s"  denotes  the  equivalent  simulation  problem.  The  ratio 
Eg/E0  can  be  made  constant  by  requiring  that 


ns  “  n 

ksas  «  k0a  (14) 

ksbs  s  kob 


end  consequently  defines  the  scaled  plasma  parameters. 
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The  results  of  equation  (14)  can  also  be  obtained  directly  from 
Maxwell's  equations  (Stratton,  194lJ.  It  suffices  to  require  that  the 
two  ratios  remain  invariant 

ci  *  \ie  (  -  ) 

T 

/ 

c2=ncr  - 

where  t  is  the  characteristic  period  and  i  denotes  length.  In  the  present 
case  we  are  dealing  with  lossless  dielectrics  hence  only  c^  need  to  be 

considered.  Now  the  invariance  of  c^  requires 

2  2  2  2 
Usesis^s  14  \iZl  * 

which  leads  directly  to  the  latter  two  expressions  of  (14) .  Moreover,  it 
is  to  be  noted  that  if  the  dielectric  permittivities  are  scaled  by  the 
same  factor,  say  c^  ,  then  this  factor  must  apply  to  all  regions  of  space. 
Thus,  if  e-^  refers  to  the  plasma  layer  and  e2  to  the  outer  region  then 

"v. 

it  follows  that 


€ls  e2s 


and  consequently 


which  is  equivalent  to  first  expression  of  (14). 

The  study  of  a  table  of  dielectric  materials  [Von  Hippel,  1958] 


at  10  Gc  has  revealed  that  the  requirements  of  the  free  space  simulating 
medium,  i.e.,  low  loss -tangent,  non-corrosiveness,  and  stability  are 
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satisfactorily  met  by  Aroclor  1232  (Monsanto  Chemical) .  The  combination 
of  air  and  Aroclor  1232  (er  =2.78  and  tan  5  =  0.008)  will  result  in  an 
ns  =0.60  corresponding  to  an  electron  density  per  cubic  centimeter 
N  =  7.83  x  10^  fQ^  where  fQ  is  the  actual  antenna  operating  frequency 
in  Gc.  With  the  simulating  frequency  fs  =  10  Gc,  air  layer  as  =  2.9  cm, 
slot  mean  radius  bs  =  0.675  cm  and  the  simulating  tank  containing  Aroclor 
1232,  it  follows  that  a/XQ  =  1.612  and  b/XQ  =  0.375.  As  a  consequence  of 
the  scaling  defined  in  (14)  and  since  the  plasma's  index  of  refraction  is 
a  function  of  the  wave  frequency,  the  plasma  environment  that  can  be  repre¬ 
sented  by  this  system  will  depend  on  the  wave  frequency  chosen.  Table  I 
shows  typical  physical  conditions  of  plasma  sheath  antenna  environment  that 
can  be  simulated  with  the  Aroclor  1232  as  free  space  and  air  as  the  plasma 
layer  dielectric. 


TABLE  I 

PLASMA  SHEATH  ANTENNA  ENVIRONMENT  REPRESENTABLE 

BY  THE  SYSTEM 


fo(Gc) 

N(elect./cm3 

a  (cm) 

b(cm) 

3 

7.02  x  1010 

16.1 

3.75 

6 

2.80  x  1011 

8.1 

1.88 

9 

6.32  x  1011 

5.4 

1.25 

12 

1.13  x  1012 

4.0 

0.94 

15 

1.75  x  1012 

3.23 

0.75 

-  ft 
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A  semicylindrical  tank  of  22  inches  inside  radius  and  24  inches 
high,  as  shown  in  Figure  3,  was  designed  and  subsequently  built  to  the 

design  specifications  by  the  Paramount  Plastic  Fabricators  of  Downey,  ; 

,  1 

California.  The  wall  material  is  plexiglass,  er  «  2.59.  The  flat  wall  of 

the  tank  is  made  of  1  inch  plexiglass  plate  and  the  curved  wall  of  0.25 

inch  plexiglass.  The  tank  holds  78  gallons  of  Aroclor  1232  oil?er  *  2.78.  ; 

j 

At  the  operating  frequency  of  10  Gc,  the  tank  allows  a  separation  of  the 

transmitting  and  the  receiving  antennas  by  at  least  28  wavelengths.  ■ 

Since  the  relative  dielectric  constants  of  the  tank  material  and  j 

l 

the  oil  are  so  close  in  value,  the  interface  between  them  has  negligible 
effect  on  the  wave  propagation.  The  reflections  from  the  curved  plexiglass- 
air  interface  were  successfully  reduced  to  a  desirable  level  by  placing  high 
performance  microwave  absorber  against  the  outside  wall  of  the  tank. 

A  ground  plane  was  constructed  using  0.312-inch  thick  aluminum 
plate  mounted  on  0.75-inch  thick  plywood.  A  2-inch  diameter  hole  was  bored 
in  the  center  of  the  plate  and  then  counter-bored  to  2.5-inch  diameter.  The 
counter-bored  opening  was  accurately  fitted  with  a  separate  machined  piece 
which  contained  the  annular  slot  antenna. 
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With  the  simulation  tank  located  on  top  of  a  stationary  bench, 
the  ground  plane  could  be  situated  conveniently  at  arty  distance  in  back 
of  the  simulation  tank  thus  allowing  the  simulation  of  the  plasma  sheath 
of  arbitrary  thickness.  The  ground  plane  was  mounted  in  the  vertical  plane 
atop  a  wooden  cart  equipped  with  rollers  for  ease  of  movement.  The  gap 
between  the  ground  plane  and  the  simulation  tank  was  adjusted  to  the  desired 
dimension  by  means  of  spacers  and  clamps  to  avoid  accidental  movement.  The 
effects  of  the  finite  dimensions  of  both  the  ground  plane  and  the  simulated 
plasma  layer  were  further  aleviated  by  placing  high  performance  microwave 
absorber  around  the  finite  boundaries  of  the  simulated  plasma  layer.  The 
absorber  in  this  position  acts  as  a  matched  termination  of  the  plasma  layer 
edges . 

In  the  early  simulation  attempt,  the  ground  plane  was  situated 
inside  the  simulation  tank  and  a  sealed  styrofoam  layer  was  used  for  the 
simulation  ox  the  plasma.  However,  this  arrangement  proved  inadequate  because 
both  the  styrofoam  and  the  ground  plane  had  to  be  rather  small  in  size  thus 
producing  considerable  open-ended  waveguide  radiation.  Difficulties  were 
also  encountered  in  the  sealing  of  the  styrofoam  sandwich. 

The  pick-up  antenna  consisted  of  a  half-wave  electric  dipole  on 
the  end  of  a  RG58A/U  coaxial  cable.  The  definition  of  the  plane  of  the  dipole 
was  achieved  by  embedding  the  dipole  in  a  plexiglass  strip  while  its  feeding 
cable  was  taken  out  of  the  tank  along  an  L~bracket  also  made  of  plexiglass. 

The  L-bracket  was  rigidly  connected  to  a  vertical  shaft,  the  axis  of  which 
coincided  with  the  center  of  the  aperture  in  the  ground  plane.  This  type  of 
arrangement  allowed  approximately  +70°  coverage  in  the  \Jr  direction  for  the 
largest  simulated  plasma  layer  of  2-inches  in  accordance  with  the  slot  co¬ 
ordinates  as  shown  in  Figure  1. 
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The  output  of  the  signal  source  was  optimized  at  the  operating 
frequency  to  maximize  the  sensitivity  of  the  system  and  frequency  stabili¬ 
zation  was  maintained  during  the  measurements.  The  signal  generator  power 
was  supplied  by  the  Varian  X-13  tube  at  10  Gc.  Frequency  stabilization 
was  achieved  through  the  Dymec  2650A  oscillator  synchronizer.  A  Varian 
V-262  klystron  was  utilized  as  the  local  oscillator  and  it  was  operated  at 
9.970  Gc.  Submersion  of  the  tube  in  an  oil  bath  assured  the  required 
frequency  stability.  The  local  oscillator  was  modulated  with  1000  cps 
square  wave.  The  signals  from  the  pick-up  dipole  and  the  local  oscillator 
were  fed  into  the  balanced  mixer-preamplifier  unit  producing  a  30  me  output 
which  was  subsequently  fed  into  the  high  gain  IF  amplifier  for  further  amp¬ 
lification.  The  1000  cps  square  wave  provided  satisfactory  modulation  of 
the  microwave  signal  so  that  an  audio  amplifier  could  be  used  with  a  crystal 
or  a  bolometer  for  detection,  aiding  in  the  calibration  of  the  receiver.  The 
receiver  output  was  fed  into  the  Antlab  1966  antenna  pattern  recorder  which 
is  equipped  with  a  selective  filter  amplifier  with  1000  cps  center  frequency 
and  4,  20,  and  50  cycle  bandwidths  which  feeds  the  pen  system  having  linear, 
logarithmic,  and  square  root  modes . 
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The  availability  of  theoretical  radiation  patterns  of  an  annular 
slot  [lyras.  1962]  and  a  thin,  long  rectangular  slot  [Omura,  1962]  motivated 
the  experimental  investigation  of  these  two  types  of  radiators.  An  annular 
slot  was  constructed  using  brass  stock  because  of  the  ease  for  precision 
machining.  Since  an  annular  slot  of  such  narrow  width  has  a  very  small  im¬ 
pedance,  a  matching  transformer  was  necessary  to  match  the  slot  impedance 
to  the  50-ohm  coaxial  feed  line.  A  coaxial  conical  taper  of  5.2  wavelengths 
was  machined  as  an  integral  part  of  the  slot  piece  and  fulfilled  the  imped¬ 
ance-matching  requirements.  Nylon  with  er  *  3.02  at  10  Gc  served  as  the  di¬ 
electric  of  the  coaxial  conical  taper.  The  desired  slot  has  a  mean  radius 
of  0.256-inch  and  0.020-inch  width.  Slots  of  smaller  width  are  practical 
but  their  impedance  tends  to  zero  and  minimizes  their  efficiency  as  radiators. 
Figure  4  shows  the  schematic  of  the  slot  and  tapering  section  with  detailed 
machining  dimensions.  The  rectangular  slot  was  constructed  using  two  wave¬ 
guide  sections  as  shown  in  Figure  5.  Two  brass  pieces  were  machined  and  in¬ 
serted  in  the  guides,  with  the  small  dimension  of  the  guide  gradually  taper- 
in  to  a  0.0312-inch  slot. 
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Because  of  the  inherent  difficulties  encountered  in  the  measurement 
of  the  slot  antenna  admittance  at  the  aperture  itself,  the  VSWR  and  the  shift 
of  the  null  from  the  short  circuit  to  the  loaded  condition  were  measured  at 
the  input  of  the  tapered  transition  section.  These  two  parameters  were  mea¬ 
sured  by  means  of  a  standard  slotted-section  technique  for  different  simulated 
plasma  thicknesses.  The  knowledge  of  the  VSWR  and  the  null  shift  permits  the 
calculation  of  the  slot  admittance  by  using  the  Smith  Chart  or  the  well-known 
transmission  line  equation.  The  measured  value  of  the  admittance  could  be 
transferred  to  any  point  on  the  taper  had  the  variation  of  the  characteristic 
admittance  with  the  taper  length  been  known  which  is  not  the  case  for  the 
linear  taper  used.  A  short-circuit  test  performed  on  the  taper  indicates  that 
this  structure  has  negligible  losses  and,  consequently,  if  can  be  considered 
as  an  ideal  transformer.  Hence,  the  functional  dependence  of  the  layer  thick¬ 
ness  on  the  slot  radiation  admittance  as  measured  at  the  input  cf  the  taper 
should  be  the  same  as  it  would  be  at  the  plane  of  the  slot  itself. 

3.2  Results 

The  experimental  antennas  patterns  are  shown  in  Fig.  6  for  the  annular 
slot.  Of  importance  to  note  is  the  occurrence  of  the  radiation  maxima  near  the 
theoretical  critical  angle,  9C  *=  37°  and  a  low  gain  region  near  9  =*  0°  .  This 
is  entirely  as  predicted  by  the  theory  [Tyras,  1962]  .  Theoretical  curves 
for  a  «  .353  XQ  and  for  a  =  1.41  \Q  have  been  plotted  and  are  compared  with 
the  experimental  values  in  Figs.  6a  and  6b  respectively.  Figure  6c  shows 
experimental  curves  for  a  =*  .706  \0  and  for  a  =  1.06  \0.  Very  good  agreement 
with  theory  is  evident  in  Fig.  6a,  but  significant  discrepancies  are  noted 
in  Fig.  6b.  As  will  be  seen  later,  the  same  situation  occurred  in  the  case 
of  the  rectangular  slot. 
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The  results  of  the  admittance  measurements  for  the  annular  slot 
are  seen  in  Fig.  7.  Reference  to  theoretical  curves  [Galejs,  19 64]}  shows 
that  the  experimental  curves  are  of  the  same  general  shapes  as  those  predict¬ 
ed.  One  notable  discrepancy  is  seen  in  the  curve  for  the  normalized  conduc¬ 
tance.  The  theory  predicts  approximately  a  5:1  ratio  between  the  maximum 
and  minimum  values  of  G/Y0,  while  the  experimental  curve  shows  approximately 
a  1.1:1  ratio. 

The  radiation  patterns  of  the  rectangular  slot  are  seen  in  Fig.  8. 
Again,  note  the  occurrence  of  radiation  maxima  near  the  predicted  critical 
angle,  0C  *  37°  .  Theoretical  curves  for  a  =  .353  X  and  for  a  *  1.41  X0 
have  been  plotted  and  compared  with  the  experimental  values  in  Figs.  8a  and 
8b  respectively.  As  observed  previously,  the  pattern  for  the  small  sheath 
thickness  of  .353  X0  agrees  quite  well  with  theory,  while  the  same  signifi¬ 
cant  discrepancies  are  noted  for  the  large  thickness.  Figure  8c  shows  ex¬ 
perimental  radiation  patterns  for  sheath  thicknesses  of  .706  XD  and  1.06  X0. 

In  evaluating  these  results,  it  is  well  to  keep  in  mind  any  approxi¬ 
mations  made  in  obtaining  the  theoretical  radiation  patterns.  For  the  problem 
at  hand,  the  method  of  saddle  point  integration  has  been  used.  This  method 
allows  the  representation  of  an  integral  in  the  form  of  a  series  of  inverse 
powers  of  distance  from  the  source  [Brekhovskikh,  I960]  .  It  is  assumed  that 
a  portion  of  the  integrand,  which  in  the  present  problem  is  the  angular  vari¬ 
ation  of  the  pattern,  varies  slowly.  That  is,  its  derivatives  are  sufficiently 
small  so  that  the  expansion  may  be  limited  to  the  first  terra.  This  approxima¬ 
tion  becomes  better  further  from  the  source. 

Now,  reference  to  Fig.  8  shows  that  for  large  plasma  thicknesses, 
the  angular  variation  becomes  quite  rapid  near  the  critical  angle.  In  fact, 
one  can  show  that  at  the  critical  angle,  the  first  derivative  is  proportional 
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to  the  square  of  the  layer  thickness.  The  situation  then,  is  the  following. 

In  the  case  of  small  plasma  thicknesses,  the  slowly  varying  assumption  is 
valid  for  the  radial  spacing  between  source  and  pickup  available  to  the  simu¬ 
lation  technique.  As  the  thickness  increases,  however,  the  angular  variation 
near  the  critical  angle  becomes  more  rapid  and  the  calculated  results  become 
less  accurate.  Physical  limitations  prevent  the  increase  of  radial  spacing. 
Thus,  to  obtain  better  agreement  with  the  simulated  patterns,  it  becomes 
necessary  to  include  a  second  term  in  the  saddle  point  approximation.  At  pre¬ 
sent,  however,  the  second  term  has  not  yet  been  evaluated  and,  consequently, 
the  theoretical  results  include  the  contribution  of  the  first  term  only. 

Another  major  factor  contribution  to  the  discrepancies  in  the  radi¬ 
ation  patterns  stems  from  the  fact  that  the  physical  boundaries  of  the  system 
are  not  sufficient  to  provide  the  minimum  required  far-field  distance  for 
thicknesses  greater  than  0.45  XQ.  The  dependence  of  the  minimum  far-field 
distance  on  the  sheath  thickness  is  seen  in  Fig.  9.  This  curve  was  plotted 
by  considering  the  half  power  bandwidth  of  the  radiation  pattern  around  the 
critical  angle.  The  half  power  bandwidths  were  calculated  from  equation  (lib). 
Now  from  Silver  [l94S]  we  have 

0HPBW  *  60°  ^  min  > 

which  yields  d/\  ,  and  subsequently  the  minimum  far-field  distance.  Thus 
y  the  radiation  pattern  for  a  <*  .353  \0  is  well  in  the  radiation  field,  while 

that  for  a  «  1.41  X0  is  in  the  near  field. 

The  admittance  curves  for  the  rectangular  slot  are  seen  in  Fig.  10. 
Theoretical  calculations  have  been  made  by  Galejs  [l964].  Comparison  with 
these  curves  shows  that  the  experimental  curve  for  the  normalized  conductance 
agrees  quite  well  with  the  theory.  The  same  can  be  said  for  the  normalized 
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susceptance  curve  except  for  very  thin  layers. 

4.  Conclusions 

In  the  present  work  only  the  simplest  case  of  the  plasma  sheath 
environment  was  treated,  i.e.,  the  case  of  a  lossless,  isotropic  and  homo¬ 
geneous  plasma.  The  simulation  technique  has  been  successfully  tested  only 
for  this  simplest  case  but  it  is  clear  that  it  may  be  used  for  the  study  of 
more  complex  plasma  sheath  models.  The  inherent  simplicity  of  this  experi¬ 
mental  system  and  its  versatility  suggest  that  the  method  is  capable  of 
simulating  more  realistic  plasma  sheath  environments  which  are  not  readily 
amenable  to  analytic  solutions. 

The  most  serious  limitation  of  the  present  technique  is  the  fact 
that  the  method  is  not  capable  of  simulating  plasmas  characterized  by  a 
negative  dielectric  constant.  Moreover,  the  effective  simulated  dielectric 
constant  of  0.36  achieved  with  the  Aroclor  1232  -  air  combination  cannot  be 
readily  lowered  because  of  the  lack  of  suitable  liquid  dielectrics  with 
>  2.8.  However,  since  the  reentry  communication  systems  can  be  expected 
to  work  well  above  the  plasma  frequency,  the  above  limitations  may  not  be 
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Fig.  6a.  Experimental  and  Theoretical  Radiation  Patterns  of 
the  Plasma  Clad  Annular  Slot 
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